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Abstract 

It is shown that it is possible to account for all three experimental 
indications for neutrino oscillations with just three neutrino flavors. In 
particular, we suggest that the solar and atmospheric neutrino anoma- 
lies are to be explained by the same mass difference and mixing. Possi- 
ble implications and future tests of the resulting mass-mixing pattern 
are given. 

1 Introduction 

Currently there are three pieces of evidence which suggest that neutrinos 
have non-zero mass differences and mixings. These are: (i) the observations 
of solar neutrinos, (ii) the anomaly in the v^j v e ratio in atmospheric neutrinos 
at low energies and (iii) the possible — V e conversion seen in the LSND 
experiment. With the conventional interpretation of these effects as being 
due to neutrino oscillations; the solar neutrino anomaly needs a 5m 2 {v e — v x ) 
of either about lCr 5 -lCr 6 eV 2 (MSW) or about lCT 10 eV 2 (long wavelength 
vacuum oscillations), the atmospheric neutrino anomaly calls for a 5m 2 (v e — 
Vp) or (ujj, — v T ) of 1CT 2 — 10~ 3 eV 2 and the LSND effect needs a 5m 2 [y e — v^) 
in the neighborhood of 1 — 2 eV 2 . For these three independent Sm 2 's at least 
one more neutrino state (beyond the three flavors) is necessary Q. 

In this letter we explore the possibility that all the neutrino anomalies 
may yet be accounted for with just three flavors of neutrinos. We assume 
only two distinct values of 5m 2 's. One value of 5m 2 is selected to explain 
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the low energy atmospheric data, while the second value is selected with the 
LSND effect in mind. 

Specifically we choose the following spectrum of Sm 2 's: 

5m 2 31 ~ 5m 2 32 ~ (1 - 2)eV 2 (1) 

5m 2 21 ~ 1(T 2 eF 2 (2) 

We then seek to determine if, with this spectrum of Sm 2 's, an explanation 
of the LSND, solar, atmospheric neutrino data can be found by appropriate 
choice of neutrino mixing angles. 

We begin by calculating the neutrino survival and transition probabilities. 
In general, these are given by 



a/3 



J2U pi eM-iEit)Ki\ 2 (3) 



Here the U a i are elements of the matrix U describing the mixing between the 
flavor eigenstates (u a ) and the mass eigenstates (i^); that is u a = Y^iU a iUi. 
For now we ignore possible CP violation, then U is real, and Eq (|3]) may be 
written as 

Pap = E^)^") 2 + 2j2 u PiUp j u ai u aj x cos y-i^j ( 4 ) 

where 5m% = m 2 — m 2 and L is the distance between the neutrino source 
and detection. We present below an explicit form of the 3x3 matrix U 

12^13 < -^13'- , 12 <->13 

U = | — C23S12 — S23S13C12 C23C12 — S23S13S12 S23C13 j (5) 
23^12 — <^23'- , 13 ( ^12 — <J23^12 _ <^23>- , 13>- , 12 <^23^13 

where C\i = cos #12, S12 = sin #12, etc.. The explicit form of the transition 
probabilities depends on the spectrum of the <5m 2, s. For the choice of <5m 2 's 
considered here, all of the oscillating terms in Eq (|]) average to zero for 
the energies and path lengths relevant to both low energy atmospheric and 
solar neutrinos. Hence, for our model, the form of the transition and survival 
probabilities relevant to solar and atmospheric neutrinos are: 

Pee = E(^) 4 ( 6 ) 
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*w = E(^) 4 (7) 

i 

Pe» = P»e = YK U eiU»i) 2 (8) 
% 

Note that the above expressions are functions of the mixing angles only, and 
are independent of the neutrino energy. 



2 Solar Neutrinos 

The four currently operating solar neutrino experiments report the fol- 
lowing results: 



Experiment 
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Recently measurements of the reaction 7 ± 8 B — > 7 Be ± p have been 
made[Q. These suggest that the cross-section for the inverse reaction 7 Be + 
p — > 8 B + 7 at energies relevant for the solar core may be somewhat smaller 
than the value used in the Standard Solar Model(SSM) of Bahcall et al. 
calculations. Hence it is possible that the flux of 8 B neutrinos is somewhat 
smaller than the SSM, while the other neutrino fluxes are unaffected. We 
allow for this possibility by defining fs as 

fB = ^P (9) 

where is the 8 B neutrino flux predicted in the SSM of Bahcall and 
Pinsonneault, which incorporates helium diffusion 0. [] Thus the parameter 
fs describes the deviation of the actual 8 B neutrino flux from the SSM value. 

We can now proceed to describe the expected counting rates at the var- 
ious solar neutrino experiments in terms of Jb and P ee the solar u e survival 
probability. 

1 Solar models which incorporate helium diffusion yield values for the depth of the 
convective zone and primordial helium abundance which are in excellent agreement with 
helioseismological datafq]. 
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With a threshold energy of 7.5 MeV the Kamiokande water Cerenkov 
detector is sensitive only to 8 B neutrinos. The expected flux is given by: 

R(KII) = (P ee + a(l - P ee ))fB x 5.69 x lO^m-V 1 (10) 

Where 5.69 x 10 10 m _2 s _1 is the SSM prediction and a (approximately 0.16) is 
the ratio of the u^ T ) — e to v e — e scattering cross sections integrated over the 
8 B neutrino spectrum. Note that we have ignored the possibility of oscillation 
into sterile flavors and assumed that solar neutrinos not interacting as z/ e 's 
interact as either z/^'s or z/ T 's with probability 1 — P ee . 

The expected counting rate in the Homestake 37 CI experiment is given 

by 

R( 37 Cl) = (6.2/b + 1.8) x P ee SNU (11) 

where 6.2 SNU's is the expected contribution from 8 B neutrinos and 1.8 
SNU's is the contribution from all other solar neutrino fluxes. Similarly, the 
expected counting rate in the 7l Ga experiments SAGE and Gallex is given 
by 

R( 71 Ga) = (13.8/b + H7.6) x P ee SNU (12) 

where 13.8 SNU's is the expected contribution from 8 B neutrinos and 117.6 
SNU's is the contribution from all other solar neutrino fluxes. 

The results of a chi-squared analysis are shown in Fig. 1. We find that 
there is a solution at the 90% CL. when the electron neutrino survival prob- 
ability is in the range 0.4 < P ee < 0.55 and the 8 B neutrino flux is in the 
range 0.55 < fs < 0.8. This is consistent with the variation in 8 B neutrino 
flux found in an analysis of solar models f|. 

The results in Fig. 1 can be interpreted in terms of the neutrino mixing 
matrix U. From Eq( |6|) the v e survival probability P ee is a function of 612 and 
#13 only. Each allowed value of fs in Fig. 1 corresponds to an allowed range 
of P ee and hence to an allowed range of 6*i 2 and 6*i 3 . In Fig. 2 we present 
a plot of the allowed values (90% CL.) of sm(9i 2 ) and sin(6 , 13 ) for f B = 0.8 
and fs = 0.65. At f B = 0.8, P ee is required to be ~ 0.43, this can be realized 
only in three flavor mixing. Hence, as shown in Fig. 2, S12 and S13 must 

2 It was first suggested by Acker et aZ.Q that solar neutrinos might be accounted for 
by the same 5m 2 and mixing as atmospheric neutrinos and hence should show an energy 
independent suppression. A search for an energy independent fit with varying fg was 
made in Ref. |Tc[ . 
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both be nonzero. At Jb = 0.65, P ee can be greater then 0.5, which can be 
accomplished in effective two flavor mixing. As shown in the figure, there 
are allowed regions with sin(# 12 ) or sin(# 13 ) equal to zero, corresponding to 
pure v e — v T or u e — mixing respectively. 



3 LSND 

The Liquid Scintillation Neutrino Detector (LSND) experiment at Los 
Alamos reports to have observed the possible appearance of V e in an initial 
beam of z/^'s These results have been interpreted as evidence of neutrino 
oscillations and the preferred range of 5m 2 and sin 2 (2$) in a two flavor mixing 
scenario given. For definiteness, we choose sin 2 (26') ~ 1.2 x 10 -3 and 5m 2 ~ 
2 eV 2 , which lie in this range [I2|. 



In the range of L/E covered by the LSND set-up, 5mf 2 L/AE ~ and 
dm^L/AE ~ 5m\ 2 L / '4E . Then the — > V e conversion probability is given 
by 

P w = W%U 2 e , Bm*(6rr&L/4E) (13) 
Thus the three-flavor interpretation of the LSND result is obtained by letting 

sm 2 (26 LSN D) ^ MU^\ 2 \U e3 \ 2 (14) 

This may be expressed as a constraint on the three flavor mixing angles. 
Using 

\U e3 \ 2 = sm 2 (6 n ) 
\U^\ 2 = cos 2 (6 n ) S m 2 (6 23 ) (15) 



We obtain 



sin 2 (20 1 
4sin 2 (# 13 )cos 2 (# 13 ) 



Choosing Sm^ ~ Sm 2 2 to be near 2 eV 2 , the LSND results then give 
sin 2 (28 LS nd) ~ 1-2 x 10~ 3 , and we have 

• //> ^2 L2 x 10-3 

4sm 2 cos 2 0i3 
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In order to determine the range of validity of this result, constraints from 
reactor and accelerator experiments must be taken into consideration. For 
this range of 5m 2 , reactor experiments give the bound JIB[: 



U e3 \ 2 < 0.02 



and accelerator experiments give the bound[14|: 



IM 2 < °- 018 - 



(18) 



(19) 



As | C/ M 3 1 is related to \U e3 \ through Eq ( |P7D these two upper limits can be 

lave: 

(20) 



combined to form bounds on the allowed values of \U e3 \ and lU^]. We have: 



0.129 < \U e3 \ < 0.141 



and 



0.123 < \U^\ < 0.134 



(21) 



Hence the requirement that the LSND results be consistent with existing 
bounds on neutrino mixing leads to rather stringent limits on the allowed 
values of \U^\ and \U^\. We will find these constraints particularly useful 
when interpreting the atmospheric neutrino data. 



4 Atmospheric Neutrinos 

Experimentally measured atmospheric neutrino fluxes are often described in 
terms of an (observed to predicted) 'ratio of ratios' R, where 

_ (*^t/^e) observed (22) 
{y^/ z/ e)MonteCarlo 

The final results f||, [31| from Kamiokande for the low energy atmospheric 
neutrino v^/v e ratio place R at 0.62 ± 0.06 ± 0.06. The results from 1MB 
[ |j~6|j are in excellent agreement with these results. Results from non-water- 
Cerenkov detectors are somewhat varying: Soudan |I7| finds an R of 0.72 ± 
0.19 

iao7 whereas the results from Nusex |L8| and Frejus jnj are consistent 



with an R of unity, although with smaller statistics than the two large water- 
Cerenkov detectors. 
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It has been known for some time that this low energy atmospheric neu- 
trino anomaly can be explained by neutrino oscillations. For 5m 2 in the 
range (4 x 10~ 3 — 2 x 10~ 2 ) eV 2 it has been shown [^(J that this anomaly 
can be explained by — v T oscillations for 

0.6 < sin 2 (2# M _ r ) < 1.0 (23) 
or by — v e oscillations for 

0.5 < sin 2 (2^_ e ) < 1.0. (24) 
Expressing these bounds in terms of the U a i we have: 

0.3 < (P MT = J2( U ^ U n) 2 ) < 0-5 (25) 

i 

for Up — v T oscillations and 

0.25 < (P, e = ^(U^iUei?) < 0.5 (26) 

i 

for Up — u e oscillations. 

In the narrow range of \U e s\ values permitted by LSND, reactor and ac- 
celerator data (Eq |2(J) P^ T is less then 0.05 and thus inconsistent with Eq 
(p5|) ; while P^ e can take on values up to 0.48. Hence, in this region, the at- 
mospheric neutrino anomaly must be explained almost exclusively by — v e 
mixing. 

With ^23 constrained in terms of 6*13 by Eq ([T7|) and sin(#i 3 ) bound by 
Eq (|2~0D , we find that 0.25 < P^ e < 0.5 if sin(0 12 ) is in the range: 

0.38 < sin(0 12 ) < 0.92 (27) 

Thus we can express this explanation of the low energy atmospheric neutrino 
anomaly consistent with the LSND, reactor, accelerator data as the region 
of the sin(6 , 12 ) — sin(6 , 13 ) plane bounded by Eq QT7D and Eq @D). 

5 A Combined Solution to the Solar, Atmo- 
spheric and LSND Data 

It is now a straightforward matter to identify simultaneous solutions to the 
Solar, Atmospheric and LSND neutrino data. As the energy independent 
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solution to the solar neutrino problem, shown in Fig. 2, and the combined 
LSND and atmospheric neutrino solution are both expressed as regions of 
the sin(6 , 12 ) — sin(6 ) i 3 ) plane, any intersection between the allowed regions 
represents the desired solution. 

Fig. 3 presents a plot of the intersecting regions of the Solar neutrino and 
Atmospheric-LSND solutions. Fig. 3a assumes the 8 B solar neutrino flux, 
Jb-, is at 80% of its SSM value, Fig. 3b assumes fs is at 70% of its SSM value 
and Fig. 3c assumes /# is at 65% of its SSM value. There is no intersection 
in Fig. 3a and narrow region of overlap in Fig. 3b broadening somewhat in 
Fig. 3c as the 8 B neutrino suppression is allowed to increase to 0.65. 

It should be noted that the selection of any region of the sin(0i2) — sin(#i3) 
plane determines the complete set of mixing angles, and hence the neutrino 
mixing matrix U, as sin(# 23 ) is fixed by Eq (|1?D. Specifically the intersection 
region of Fig. 3b corresponds to sin(# 12 ) ~ 0.707, sin(# 13 ) ~ 0.140 and sin(# 23 ) 
~ 0.125. 

Using Eq (H) we present below the explicit form of the 3x3 mixing matrix 
U corresponding to solution region of Fig. 3b: 



While for the solution region corresponding to Fig. 3c we find that, in addi- 
tion to Eq (|28|) above, the following range of matrix values are allowed: 

/ .630 .764 .140 \ / .764 .630 .140 \ 

-.776 .619 .124 <-> -.645 .754 .124 . (29) 
\ -.010 -.187 .982 / V -.028 -.185 .982 / 



6 Implications 

(i) Both Super-Kamiokande pi] and SNO should see NO spectrum 
distortion in either v — e scattering or the v e D charged current. The suppres- 
sion in v— e scattering should be in the range 0.38 - 0.40 of SSM at all energies 
and in v e D, charged current suppression should be about fsPee ~ 0.32 — 0.34. 




(28) 
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(ii) Borexino [23| should observe the 7 Be line at a rate of [P ee + — P ee )] 



where (3 is the ratio of v^ T ) — e to v e — e scattering cross sections. We thus 
expect 0.56 to 0.58 of the SSM rate, and an identical suppression should hold 
for the pep line. 

(iii) The atmospheric v^/v e anomaly should be confirmed by Super- 
Kamiokande. Zenith angle dependence of multi-GeV neutrinos should con- 
firm the tentative evidence seen in Kamiokande |2(J . But most important is 
our prediction || that — u e oscillations should be confirmed by observation 
of excess high energy e-like upcoming shower events (above and beyond 
neutral current events). 



(iv) Future reactor experiments such as CHOOZ |24| and Palo Verde 
[ P5fl which will be sensitive to 5m 2 upto 10 -3 eV 2 should see a V e survival 
probability of P ee ~ 0.48 - 0.5. 

(v) Long baseline experiments (such as MINOS |Z6J, CERN-LNGS g7| 
and KEK-PS E362 p|) which will probe 5m 2 upto 10" 3 eV 2 should see 



— v T conversion with P^ T = Y,i(Ufj,iU T i) ~ 0.028 — 0.035 accompanied by 

— v e conversion at P^ e ~ 0.46 — 0.48. 

(vi) Short baseline experiments such as CHORUS Q, NOMAD Q and 
COSMOS |3l] will probe — v T conversion for Sm 2 > 0.1 eV 2 . We predict, 
at Sm 2 = 1 - 2 eV 2 , an effective sin 2 (2#) of 4(^ t3 [/ T3 ) 2 which is 0.06. 

(vii) Large mixings in some ranges of 6m 2 lead to strong conversion of 
to V T due to the MSW effect in supernova, leading to a harder energy 

spectrum of the emerging 77 e 's. This can lead to potential conflict with ob- 
servation of neutrinos from SN1987A. For the 5m 2 in our scenario this is not 



a problem |32 



(viii) The neutrino mass spectrum implied by our scenario is: 

rrii ~ m o 
ni2 ~ mo + e 



m 3 ~ sjml + 2eV 2 (30) 

where 5m 2 2 ~ (2moe+e 2 ) ~ 10~ 2 eV 2 . There are two limiting cases of interest 
assuming that the largest mass is in the eV range. One is the hierarchical 
limit, in which mo is negligible. Then m-i m 2 ~ 0.1e\^ and m 3 ~ lAeV. 
The other is the nearly degenerate limit, in which 

m-i ~ leV 
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m 2 ~ (1 + e)eV 

m 3 ~ 1.73eV (31) 

with e ~ |(10 _2 )e^. Then, the sum of the neutrino masses is J2i m i ~ AeV. 
In this case, the Cosmological density parameter associated with neutrinos 
Vt v = 0.011/i- 2 J2im = 0.044/i" 2 « 0.2 (for h of about 0.5) and the amount 
of neutrino dark matter component along with cold dark matter makes for a 
viable and testable scenario for mixed dark matter |33 . 



(ix) When the neutrinos are Majorana particles, the effective mass 
< m Ue > relevant in neutrino-less double /3-decay analysis is 

<m„ e (32) 

i 

We find that in the case of the hierarchical spectrum < m Ue >~ O.leV^ 
whereas in the degenerate case < m Ve >~ leV (this could be somewhat 
smaller when CP phases are taken into account). It is interesting that these 
values are in the range of what the double beta decay experiments can probe 
now and in the near future [33]. 

(x) When the mixing matrix is allowed to have a CP violating phase, the 
CP violating neutrino flavor conversion probability differences are given by 

m ' 



-ifiT f*T -'■fie -'fie 

= -U v cp [sin D 12 + sin D 23 + sin D 31 ] (33) 

where 

r CP = im [u^u; 2 u, 3 u: 3 } 

= IMI^IIMI^slsin^, (34) 

and 

Dij = 5m 2 lj L/2E (35) 

with being the phase in the mixing matrix. With the matrix of Eq.( 
|2TS| ), Jq P < 0.07; and [sin_D 12 + sin_D 2 3 + sin_D 31 ] »s sin_D 12 is given by ~ 
-1 for L/E = 730km/10GeV (relevant for MINOS) and also for L/E = 
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250km/3GeV (relevant for E362). Hence, AP can be as large as 0.07. (For 
these parameters, matter effects are negligible [j35|). 



We conclude by stressing that our proposal to account for both solar 
and atmospheric neutrino anomalies by the same mass and mixing can be 
confirmed or ruled out in the very near future. 
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Figure Captions 



Figure 1 Contour plot showing the allowed values of the parameters P ee 
and /en at the 90 % (solid line) and 95 % (dashed line) confidence levels in 
the three flavor mixing solution to the solar neutrino problem. 

Figure 2 Contour plot showing the allowed values of sin(#i2) and sin(#i3) 
(90 % confidence level) in the three flavor mixing solution to the solar neu- 
trino problem for fixed values of Solid line; f B = 0.8, dashed line; 
f B = 0.65. 

Figure 3 Combined solution: solar atmospheric and LSND results. Dashed 
line; bounds from LSND, reactor and accelerator experiments, Hatched re- 
gion; Atmospheric neutrino anomaly explained by — v e oscillations, Solid 
line; allowed region at the 90 % C.L. in the three flavor mixing solution to 
the solar neutrino problem for: (a) f B = 0.8 (b) f B = 0.7 and (c) f B = 0.65. 
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